Astronomy & Astrophysics manuscript no. article 


©ESO 2012 


November 6, 2012 





WASP-52b, WASP-58b, WASP-59b, and WASP-60b: 
four new transiting close-in giant planets 

G. Hebrard 1 ' 2 , A. Collier Cameron 3 , D. J. A. Brown 3 , R. F. Diaz 41 - 2 , F. Faedi 13 , B. Smalley 5 , D. R. Anderson 5 , 

D. Armstrong 13 , S. C. C. Barros 4 , J. Bento 12 , F. Bouchy 1 - 2 , A. P. Doyle 5 , B. Enoch 3 , Y. Gomez Maqueo Chew 1314 , 

E. M. Hebrard 2 , C. Hellier 5 , M. Lendl 6 , T. A. Lister 11 , R F L. Maxted 5 , J. McCormac 8 , C. Moutou 4 , D. Pollacco 13 , 

D. Queloz 6 , A. Santerne 42 , 1. Skillen 8 , J. Southworth 5 , J. Tregloan-Reed 5 , A. H. M. J. Triaud 6 , S. Udry 6 , 
M. Vanhuysse 10 , C. A. Watson 7 , R. G. West 9 , P. J. Wheatley 13 

1 Institut d' Astrophysique de Paris, UMR7095 CNRS, Universite Pierre & Marie Curie, 98bis boulevard Arago, 75014 Paris, France 
e-mail: hebrard@iap.fr 

2 Observatoire de Haute-Provence, CNRS/OAMP, 04870 Saint-Michel-l'Observatoire, France 

3 School of Physics and Astronomy, University of St Andrews, St Andrews, Fife KY16 9SS, UK 

4 Laboratoire d' Astrophysique de Marseille, Univ. de Provence, CNRS (UMR61 10), 38 r. F. Joliot Curie, 13388 Marseille cedex 13, France 

5 Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK 

6 Observatoire de Geneve, Universite de Geneve, 51 Chemin des Maillettes, 1290 Sauverny, Switzerland 

7 Astrophysics Research Centre, School of Mathematics and Physics, Queens University Belfast, University Road, Belfast BT7 INN, UK 

8 Isaac Newton Group of Telescopes, Apartado de Correos 321, 38700 Santa Cruz de Palma, Spain 

9 Department of Physics and Astronomy, University of Leicester, Leicester, LEI 7RH, UK 

10 OverSky, 47, allee des Palanques, BP 12, 33127 Saint-Jean dTllac, France 

11 Las Cumbres Observatory, 6740 Cortona Drive Suite 102, Goleta, CA 931 17, USA 

12 The Open University, Walton Hall, Milton Keynes, MK7 6AA, UK 

13 Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK 

14 Department of Physics and Astronomy, Vanderbilt University, 6301 Stevenson Center, Nashville, TN 37235, USA 
Received TBC; accepted TBC 

ABSTRACT 

We present the discovery of four new transiting hot jupiters, detected mainly from SuperWASP-North and SOPHIE observations. 
These new planets, WASP-52b, WASP-58b, WASP-59b, and WASP-60b, have orbital periods ranging from 1.7 to 7.9 days, masses 
between 0.46 and 0.94 M Jup , and radii between 0.73 and 1.49 Rj up . Their Gl to K5 dwarf host stars have V magnitudes in the range 
11.7 - 13.0. The depths of the transits are between 0.6 and 2.7%, depending on the target. With their large radii, WASP-52b and 
58b are new cases of low-density, inflated planets, whereas WASP-59b is likely to have a large, dense core. WASP-60 shows shallow 
transits. In the case of WASP-52 we also detected the Rossiter-McLaughlin anomaly via time-resolved spectroscopy of a transit. We 
measured the sky-projected obliquity A = 24°+", indicating that WASP-52b orbits in the same direction as its host star is rotating and 
that this prograde orbit is slightly misaligned with the stellar equator. These four new planetary systems increase our statistics on hot 
jupiters, and provide new targets for follow-up studies. 

Key words. Planetary systems - Techniques: radial velocities - Techniques: photometric - Stars: individual: WASP-52, WASP-58, 
WASP-59, WASP-60 



1. Introduction 

About 200 exoplanets are known today as transiting in front of 
their host stars as seen from the Earth. That population is par- 
ticularly interesting since it allows numerous studies, including 
accurate radius, mass, and density measurements, atmospheric 
studies in absorption through transits and in emission through 
occultations, dynamic analyses through possible timing varia- 
tions, and obliquity measurements. Increasing the size of that 
sample is essential in order both to improve our statistical knowl- 
edge on exoplanets, and to discover individual cases particu- 
larly well-adapted to follow-up studies. Several ground-based 
photometric programs are surveying large fields with that goal, 
including SuperWASP (Pollacco et al. |2006T > and HAT (Bakos 
et al. 2007). More recently, the dedicated space-based missions 
CoRoT (Baglin et al. |2009l > and Kepler (Borucki et al. [20T0l 
also joined that effort and have discovered new transiting exo- 



planets, in particular with longer periods and smaller radii than 
are easily detected by ground-based surveys. Because they sur- 
vey smaller fields, however, CoRoT and Kepler detect planets 
transiting mainly in front of fainter stars, typically in the mag- 
nitude range 13 < V < 16. Such faint stars make follow-up 
observations difficult. By comparison, the planets detected by 
SuperWASP or HAT transit brighter stars, typically in the range 
10 < V < 13, allowing easier and more accurate complemen- 
tary studies. Thus, ground-based surveys for detection of transit- 
ing exoplanets remain pertinent, and are complementary to their 
space-based counterparts. 

We present here the discovery of four new transiting plan- 
ets, detected with SuperWASP-North and SOPHIE. These two 
instruments allowed the detection of the first three WASP plan- 
ets (Collier Cameron et al. I2007al Pollacco et al. 2008 ), fol- 
lowed by about twenty others, including inflated hot jupiters 
such as WASP-12b (Hebb et al.|2009 ) and WASP-48b (Enoch 
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Fig. 1. Photometry and radial velocities of WASP-52. The two Fig. 2. Photometry and radial velocities of WASP-58. The two 

upper panels show the transit light curves and their residu- upper panels show the SuperWASP-N light curve and its resid- 

als from the fit. Different photometric observations are ranked uals from the fit. The three lower plots show the SOPHIE radial 

chronologically from top to bottom, with the same order as in velocities and their residuals from the Keplerian fit. 
Table [2] The three lower plots respectively show the phase- 
folded radial velocities, and their O-C residuals from the 
Keplerian fit as a function of transit phase and Julian date. 
SOPHIE HE1, HE2 and HR, and CORALIE data (see Table 

and Sect. O are plotted in red circle, blue diamonds, black et aL |2QTll, an d low-mass giant planets, e.g. WASP- 13b (Skillen 

downward triangles, and green upward triangles, respectively. et a i. |2QQ9l > and WASP-21b (Bouchy et al. [20T0T >. 

The four new planets presented here are WASP-52b, WASP- 
58b, WASP-59b, and WASP-60b. They are giant, close-in plan- 
ets, with masses ranging from 0.46 to 0.94 Mj up , orbital periods 
between 1.7 and 7.9 days, and transit depths between 0.6 and 
et al. 1201 11 1, compact sub-hot-jupiters like WASP- 1 lb (West et 2.7 %. We describe in Sect.[2]ihe observations that allowed their 
al. 2009), hot jupiters on eccentric orbits like WASP-38b (Barros discovery, their analyses in Sect. [3] and conclude in Sect. [4] 
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Fig. 3. Photometry and radial velocities of WASP-59. The two 
upper panels show the transit light curves and their residu- 
als from the fit. Different photometric observations are ranked 
chronologically from top to bottom, with the same order as in 
Table [2] The three lower plots show the SOPHIE radial veloci- 
ties and their residuals from the Keplerian fit. 



Fig. 4. Photometry and radial velocities of WASP-60. The two 
upper panels show the SuperWASP-N (in black) and NITES (in 
red) transit light curves and their residuals from the fit. The three 
lower plots show the SOPHIE radial velocities and their resid- 
uals from the Keplerian fit. SOPHIE HE1 and H2 (see Table 0] 
and Sect. 12.21 ) are plotted in red circles and blue diamonds, re- 
spectively. 



2. Observations 

2.1. Photometric detection with SuperWASP 

The coordinates and magnitudes of the stars WASP-52, 58, 59, 
and 60 are reported in the upper part of Table [TJ They were first 
identified as promising candidates host stars for transiting plan- 
ets with SuperWASP-North. Located on La Palma in the Canary 



Islands, Spain, it consists of eight Canon 200 mm //1.8 focal 
lenses coupled to e2v 2048 x 2048 pixel CCDs with 13.7" pixels 
and a field-of-view of 7.8° x 7.8°, associated with a custom-built 
photometric reduction pipeline (Pollacco et al. 2006 ). It observes 
with a broad-band custom filter (400-700 nm). Thousands of 
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Table 1. Stellar parameters of the four planet host stars from analysis of SOPHIE spectra. 



Parameter 


WASP-52 


WASP-58 


WASP-59 


WASP-60 


RA (J2000) 


23:13:58.76 


18:18:48.25 


23:18:29.54 


23:46:39.98 


DEC (J2000) 


+08:45:40.6 


+45:10:19.1 


+24:53:21.4 


+31:09:21.4 


B 


12.9 


12.33 


13.92 


12.86 


V 


12.0 


11.66 


13.00 


12.18 


7W (K) 


5000 ± 100 


5800 ± 150 


4650 ± 150 


5900 ± 100 




4.5 ±0.1 


4.3 ±0.1 


4.55 ±0.15 


4.2 ±0.1 


£(kms-') 


0.9 ±0.1 


1.1 ±0.2 


0.7 ±0.3 


1.0 ±0.1 


Vmac (kms~') 


0.7 ± 0.3 


2.3 ±0.3 





2.9 ±0.3 


vsirn, (kms~') 


3.6 ± 0.9 


2.8 ±0.9 


2.3 ± 1.2 


3.4 ±0.8 


P*» (d) 


11.8 ±3.3 


21.1 ±7.6 


15.9 ±8.4 


20.6 ± 5.5 


Gyrochron. age (Gyr) 


4 +0 3 

u - M).2 


q 9+4.5 
j.z. 2 j 


5 +0 - 7 


3 6 +43 

-2.1 


log^ K 


-4.4 ± 0.2 


-4.4 ± 0.2 


-4.1 ±0.2 


-4.4 ± 0.2 


[Fe/H] 


0.03 ±0.12 


-0.45 ± 0.09 


-0.15 ±0.11 


-0.04 ± 0.09 


[Na/H] 


0.29 ± 0.08 


-0.44 ± 0.07 


- 


0.07 ±0.11 


[Mg/H] 


0.26 ±0.13 


-0.30 ± 0.07 


- 


0.02 ± 0.06 


[Si/H] 


0.25 ±0.18 


-0.34 ± 0.06 


- 


0.15 ±0.08 


[Ca/H] 


0.04 ±0.13 


-0.30 ±0.13 


0.23 ±0.18 


0.06 ±0.13 


[Sc/H] 


0.17 ±0.11 


-0.22 ± 0.07 


- 


0.18 ± 0.08 


[Ti/H] 


0.06 ±0.12 


-0.29 ± 0.07 


0.28 ±0.18 


0.00 ± 0.07 


[V/H] 


0.19 ±0.13 




0.66 ±0.19 


0.02 ±0.10 


[Cr/H] 


0.11 ±0.15 


-0.39 ±0.18 




0.05 ±0.14 


[Co/H] 


0.20 ± 0.07 




0.16 ± 0.1 1 


-0.01 ±0.15 


[Ni/H] 


0.11 ±0.13 


-0.49 ± 0.09 


-0.10 ±0.10 


0.29 ± 0.08 


log A(Li) 


< 0.3 ± 0.2 


2.32 ±0.13 


< 0.17 ±0.21 


< 0.83 ± 0.09 


Sp. Type 


K2V 


G2V 


K5V 


G1V 


Distance (pc) 


140 ± 20 


300 ± 50 pc 


125 ± 25 


400 ± 60 



photometric points were secured with SuperWASP-North over 
several seasons for each target (see Table|2]for a summary). 

Periodic signatures of possible planetary transits were iden- 
tified in these light curves using the algorithms presented by 
Collier Cameron et al. ([2006b . Phase-folded SuperWASP light 
curves are plotted in the upper panels of Figs. Q] [2] [3] and|4] for 
WASP-52, 58, 59, and 60, respectively. They show dips of the or- 
der of one percent, compatible with transiting giant planets with 
periods of 1.7, 5.0, 7.9, and 4.3 days, respectively. 

2.2. Radial velocities with the spectrograph SOPHIE 

The spectroscopic follow-up of these four candidates was done 
with SOPHIE, the spectrograph dedicated to high-precision ra- 
dial velocity measurements at the 1.93-m telescope of Haute- 
Provence Observatory, France (Bouchy et al. [2009). The first 
goal of these observations is to establish the planetary nature 
of the transiting candidates found in photometry, which in most 
cases are not due to planets but to astrophysical false positives 
such as blended stellar binaries. The second goal is to character- 
ize the secured planets by measuring in particular their masses 
and orbital eccentricities. 

SOPHIE was mainly used in High-Efficiency mode (HE) 
with a resolving power R = 40000 in order to increase the 
throughput for these faint targets. The exposure times ranged 
from 10 to 30min depending on the targets, and they were ad- 
justed as a function of the weather conditions in order to keep 
the signal-to-noise ratio as constant as possible for any given 
target. The spectra were extracted using the SOPHIE pipeline, 
and the radial velocities were measured from weighted cross- 
correlation with numerical masks characteristic of the spectral 
type of the target (Baranne et al. 1996; Pepe et al. l20021 >. We ad- 
justed the number of spectral orders used in the cross-correlation 



in order to reduce the dispersion of the measurements. Indeed, 
some spectral domains are noisy (especially in the blue part of 
the spectra), so using them degrades the accuracy of the radial- 
velocity measurement. 

The error bars on the radial velocities were computed from 
the cross-correlation function using the method presented by 
Boisse et al. (f2010). Some spectra were contaminated by moon- 
light. Following the method described in Pollacco et al. |2008 ) 
and Hebrard et al. (2008), we estimated and corrected the moon- 
light contamination by using the second SOPHIE fibre aperture, 
which is targeted on the sky whereas the first aperture points 
toward the target. This typically results in radial velocity correc- 
tions of a few ms -1 up to a few tens of ms~'. Removing these 
points does not significantly modify the orbital solutions. 

The radial velocity measurements are reported in Tableland 
are displayed in the lower panels of Figs. HHI] together with 
their Keplerian fits and the residuals. They show variations in 
phase with the SuperWASP transit ephemeris and with semi- 
amplitudes of a few tens ms -1 , implying companion masses 
slightly below a Jupiter mass. Radial velocities measured us- 
ing different stellar masks (F0, G2, or K5) produce variations 
with similar amplitudes, so there is no evidence that their vari- 
ations could be explained by blend scenarios caused by stars of 
different spectral types. Similarly, the cross-correlation function 
bisector spans show neither variations or trends as a function 
of radial velocity (Fig. [5]). This reinforces the conclusion that 
the radial-velocity variations are not caused by spectral-line pro- 
file changes attributable to blends or stellar activity. We thus 
conclude that our four targets harbour transiting giant planets, 
which we designate as WASP-52b, WASP-58b, WASP-59b, and 
WASP-60b hereafter. 

The raw residuals of the Keplerian fits to two of the four tar- 
gets (WASP-52 and 60) show an abrupt offset of ~ 100 ms -1 
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Fig. 5. Bisector span as a function of the radial velocities with 1 -<x error bars for WASP-52, 58,59 and 60 (from left to right and top 
to bottom). Color codes are the same as in Figs. Q]|4] The ranges here have the same extents in the x- and y-axes. 



around BJD = 2 455 500 which could be due to possible extra 
companions in the planetary systems. However, the fact that the 
offset is similar in amplitude and time for both targets suggests 
it could be instrumental. WASP-58 and 59 were not observed 
at that epoch. To constrain possible extra companions we re- 
observed WASP-52 the following season with SOPHIE in High- 
Resolution mode (HR), with a resolving power R = 75 000 fol- 
lowing the installation of new octagonal fiber sections in the 
spectrograph, which improved significantly the radial-velocity 
accuracy of SOPHIE (Bouchy et al. I20T21 . WASP-52 was also 
observed with the CORALIE spectrograph mounted on the 1.2- 
m Euler-Swiss telescope at La Silla (Queloz et al. 2000). These 
new measurements did not confirm any structure in the residuals 
of the one-Keplerian fit, so we conclude that the offsets detected 
during the previous season on WASP-52 and 60 were instrumen- 
tal. Such systematic effects have already been observed in the 
past with the HE mode of SOPHIE, but not with the HR mode. 
They are not well understood, but seem to be linked with temper- 
ature variations in the telescope dome. To circumvent that effect, 
we considered the SOPHIE HE data of WASP-52 and 60 secured 
before and after the zero-point change as independent datasets, 
allowing a free radial-velocity offset to be fitted between them. 
These are noted "SOPHIE HE1" and "SOPHIE HE2" in Tableg] 



2.3. Additional photometry of the transits 

Once the radial velocities established that the four candidates 
identified with photometry actually were transiting planets, we 
obtained extra transit light curves with telescopes larger than 
SuperWASP, allowing improved resolution and more precise 
time-series photometry during transit. The goal is to refine the 
determination of the parameters derived from photometry. It was 
not possible to observe extra transits of WASP-58b with our 
available telescopes, due to the orbital period of that system be- 
ing almost exactly an integer number of days (P = 5.018 days). 

We used the six telescopes below equipped with cameras to 
observe extra transit light curves of the three remaining planets: 

- JGT: the 0.94-m James Gregory Telescope located at 
University of St. Andrews, Scotland; 

- FTS: the 2-m Faulkes Telescope South located at Siding Spring 
Observatory in New South Wales, Australia; 

- EulerCam: the new camera at the 1 .2-m Euler-Swiss telescope 
at La Silla, Chile; 

- BUSCA: at the Calar Alto Observatory 2. 2-m telescope, Spain; 

- OverSky: the 0.36-m telescope at the OverSky Observatory, La 
Palma, Spain; 

- NITES: a 0.4-m telescope at La Palma, Spain. 

The dates of the transit observations made with these facili- 
ties as well as the filters used are reported in Table [2] Standard 
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aperture photometry was carried out on all these datasets, us- 
ing comparison stars available in the different fields for relative 
photometry. The extra light curves of WASP-52, 59, and 60 are 
plotted in the upper panels of Figs. Q] [3] and[4] respectively. 



2.4. Transit spectroscopy 

A transit of WASP-52b was observed in spectroscopy with 
SOPHIE. The goal was to detect the Rossiter-McLaughlin 
anomaly, which is an apparent time-dependent distortion of the 
stellar line profiles that occurs when a planet transits the disc of 
a rotating star. It allows the measurement of the sky-projected 
angle between the planetary orbital axis and the stellar rotation 
axis, usually denoted A (see, e.g., Hebrard et al. 2008 ). The spec- 
troscopic transit was observed in the HE mode of SOPHIE to 
improve the throughput. 

The target was observed continuously during a four-hour se- 
quence under good weather conditions on 2011 August 21. 19 
measurements were secured with exposure times ranging from 
600 to 1200 seconds, including nine within the transit. The re- 
maining observations obtained before and after the transit are 
mandatory for reference points. The radial velocities were ex- 
tracted as above (Sect. 12.21 1. but using slightly fewer orders for 
the cross-correlation to reduce the dispersion of the measure- 
ments. This depends on the flux color balance in the spectra, 
which depends on the airmasses of the different observations. 

These observations are plotted in Fig. [8] and are la- 
belled "SOPHIE HE-RM" in Table H The Rossiter-McLaughlin 
anomaly is detected, with a full amplitude of ~ 80ms~'. The 
shape of the anomaly indicates a prograde orbit, with an asym- 
metry that suggests misalignment (see Sect. 13.31 ). 

For different reasons, including poor weather conditions, we 
were unable to perform time-resolved spectroscopy of the tran- 
sits of the three other planets. 



Table 2. Photometric observations. 



Telescope 


Band 


Date 


WASP-52 






SuperWASP-N 




2008 Jul - 2009 Nov 


JGT 


r 


2010 Oct 24 


FTS 


z 


2011 Aug 02 


EulerCam 


r 


2011 Aug 20 


BUSCA 


u 


2011 Aug 27 


BUSCA 


g 


2011 Aug 27 


BUSCA 


r 


2011 Aug 27 


BUSCA 


z 


2011 Aug 27 


EulerCam 


r 


2011 Sep 24 


WASP-58 






SuperWASP-N 




2004 May - 2010 Aug 


WASP-59 






SuperWASP-N 




2004 May - 2006 Nov 


OverSky 


r 


2011 Oct 28 


OverSky 


r 


2011 Nov 04 


OverSky 


r 


2011 Nov 20 


JGT 


r 


2011 Nov 20 


WASP-60 






SuperWASP-N 




2004 May - 2007 Dec 


NITES 


r 


2011 Aug 25 



3. Analysis 

3. 1 . Parameters of the host stars 

Spectral analysis of the four host stars was performed with the 
co-added individual spectra obtained with SOPHIE for radial ve- 
locity measurements (Sect. l2~2l . We used here only the spectra 
without moonlight contamination, and only the HE spectra out of 
transit in the case of WASP-52. This allows high-enough signal- 
to-noise ratio to be reached, with a given spectral resolution, and 
avoid any possible contamination in the spectra. The parame- 
ters obtained from these analyses are listed in the bottom part of 
TableQ] 

We used the methods given by |Gillon et al. (2009")| The H a 
line was used to determine the effective temperature r e ff, except 
for WASP-59 where H a was too weak; for that star T e ff was mea- 
sured from the null dependence of abundance on excitation po- 
tential. The surface gravity logg, was determined from the Cai 
lines at 6 162 A and 6439 A dBruntt et al. 201 Oak along with the 
Na i D and Mg i b lines. The spectral types were established from 
T e s using the table in Gray (2008 ). 

The elemental abundances were determined from equivalent 
width measurements of several clean and unblended lines. The 
microturbulence £ t was evaluated from Fe i using the method of 
Magain (1984) The quoted error estimates include that given by 
the uncertainties in T e g, logg, and as well as the scatter due 
to measurement and atomic data uncertainties. 

The projected stellar rotation velocities v sin /„ were deter- 
mined by fitting the profiles of several unblended Fei lines. 
They agree with the measurements secured from the cross- 
correlation function following Boisse et al. (1201 01) . The values 
for macroturbulence v mac were assumed based on the tabulations 
by Gray (2008) and [Bruntt et al. (2010b~)| Macroturbulence was 
assumed to be zero for WASP-59, since for mid-K stars it is ex- 
pected to be lower than that of thermal broadening (Gray 2008 ). 
An instrumental FWHM of 0.15 + 0.01 A was measured from 
the telluric lines around 6300 A, and best fitting values of the 
v sin U were obtained. 

Lithium is detected only in WASP-58. The temperature of 
5800 K along with the derived lithium abundance log A(Li) im- 
ply a lower age limit of ~ 2 Gyr when compared with NGC 752 
(Sest ito & Randich 20051 ). There is no significant detection of 
lithium in the spectra of WASP-52, 59 and 60, with equiv- 
alent width upper limits of 4 mA, 15 mA, and 6 mA, respec- 
tively. The lack of lithium would suggest lower age limits of 
0.5 Gy and several Gyr for WASP-52 and 60, respectively 
dSestito & Randich 20051 ). For the mid-K star WASP-59, lithium 
is expected to be depleted in only a few 100 Myr. 

For WASP-52, 58, 59, and 60, the rotation rates implied by 
the v sin/. (P m = 11.8+3.3 d, 21. 1+7.6 d, 15.9±8.4d, and20.6+ 
5.5 d, respectively) gives gyrochronological ages of 0.4^ Gyr, 
3.2ffi Gyr, 0.5 ^ Gyr, and 3.6^ Gyr, respectively, using the 
|Barnes (2007)l relation. 

WASP-52 and 59 show chromospheric emission peaks in 
the Ca ii H+K lines that indicate they are active stars, whereas 
WASP-58 and 60 do not show any evident emission. We com- 
puted the log R' HK indices following the calibration of Boisse et 
al. (120101 ) 

The SuperWASP light curve of the active star WASP-52 
shows long term modulations. EulerCam data also show irreg- 
ularities. We have analyzed the SuperWASP light curve to deter- 
mine whether they are periodic, and possibly due to the combi- 
nation of magnetic activity and the rotation of the star. The value 
of v sin U derived above together with the estimated stellar radius 
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Fig. 6. Upper panel: Periodograms for the WASP data from two seasons for WASP-52. The date range (JD-245000) is given in the 
title of each panel. Horizontal lines indicate false alarm probability levels FAP = 0.1,0.01,0.001. Lower panel: Data for the two 
seasons folded on the period P = 16.38 days. 



implies a rotation period of 5 - 15 days, assuming that the rota- 
tion axis of the star is approximately aligned with the orbital axis 
of the planet. We show below in Sect. l3.3l that this is not exactly 
the case, but the misalignment remains low. We used the sine- 
wave fitting method described in Maxted et al. (1201 11 1 to calcu- 
late the periodograms shown in Fig. [6] These are calculated over 
4096 uniformly spaced frequencies from to 2.5 cycles/day. The 
false-alarm probability levels shown in these figures are calcu- 
lated using a boot-strap Monte Carlo method also described in 
Maxted et al. (1201 It . Variability due to star spots is not expected 
to be coherent on long timescales as a consequence of the finite 
lifetime of star-spots and differential rotation in the photosphere, 
and so we analyzed the two seasons of data for WASP-52 sepa- 
rately. 

The second season of data shows a highly significant peri- 
odic modulation with a period of 16.4 + 0.04 days and an ampli- 
tude of 9.6 milli-magnitudes. We used a boot-strap Monte Carlo 
method to estimate the error on this period. We also calculated 
periodograms of three stars with similar magnitudes and col- 
ors to WASP-52 observed simultaneously with the same WASP 
camera. None of these stars showed significant periodic modu- 
lation in the same frequency interval. The strongest peak in the 
first season of data is marginally significant and is consistent 
with the 1-day alias of the period detected in the second season 
of data. 

3.2. Parameters of the planetary systems 
3.2.1. Method 

For a given system, all the available data were fitted simultane- 
ously: the radial velocities with a Keplerian orbit, and the transit 
light curves with the Mandel & Agol (2002) algorithm using 



four non-linear limb-darkening coefficients from the tabulations 
of Claret (2004). The global fit was done with the latest ver- 
sion of the Markov-chain Monte-Carlo (MCMC) procedure pre- 
sented in detail by Collier Cameron et al. (2007b ) and Pollacco et 
al. (f2008 ). The free parameters in the fits are the transit epoch To, 
the orbital period P, the depth of the transit or the planet/star area 
ratio (R p /R„) 2 , the total transit duration tj, the impact parame- 
ter b, the eccentricity e of the orbit, the argument of periastron 
o>, the semi-amplitude K\ of the radial velocity variations, the 
n systemic radial velocities y n (one for each of the n radial ve- 
locity datasets). When no significant eccentricity is detected (in 
three out of the four systems), the orbits are considered as circu- 
lar in the final fit. The stellar density p, is directly derived from 
the transit light curve (through P, tj, and b values). The density, 
together with 77 e ff and [Fe/H], yield the stellar mass M, and ra- 
dius R* through empirical calibrations established for eclipsing 
binaries (Torres et al. I2010I Enoch et al. I201 11 1. From these we 
conclude that the four host targets are dwarf stars. The plane- 
tary mass M p and radius R p are derived from K\ and (R p /R*) 2 , 
which directly provides the planetary density p p . We also com- 
pute the derived parameters as the surface gravity log g* (which 
here is more accurate than that obtained from spectral analysis), 
the semi major axis a and the orbital inclination i p , as well as 
the equilibrium temperature 7p of the planet, assuming it to be 
a black body and that energy is efficiently redistributed from the 
planets day-side to its night-side. 

For comparison, the stellar mass is also derived from mass 
tracks in p»-vs.-77 e ff HR diagrams (Sozzetti et al. 2007). We com- 
pared several models for the tracks in the HR diagrams: the 
Padova isochrones of Marigo et al. (f2pl0), the Teramo stellar 
models from Pietrinferni et al. (120041 ). the Victoria-Regina stel- 
lar models from Vandenberg et al. (120061 ). and the Yonsei-Yale 
isochrones of Demarque et al. ( 120041 ) . The last among these are 
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Table 3. Planetary and stellar parameters for the systems WASP-52, 58, 59, and 60. 



Parameter 


Symbol 


WASP-52 


WASP-58 


Units 


Transit epoch (HJD-2450000.0) 


To 


5793.68143 ±0.00009 


5183.9335 ±0.0010 


days 


Orbital period 


P 


1.7497798 ±0.0000012 


5.017180 ±0.000011 


days 


Transit duration 


h 


0.0754 ± 0.0005 


0.1582 ±0.0043 


days 


Planet/star area ratio 


(RJR,) 2 


0.0271 ±0.0004 


0.0145 ±0.0010 




Impact parameter 


b 


0.60 ± 0.02 


0.46 ± 0.23 


R, 


Scaled stellar radius 


RJa 


0.1355 ±0.0020 


0.097 ±0.011 




Stellar density 


p. 


1.76 ± 0.08 


0.58 ±0.19 


Po 


Stellar surface gravity 


l°gg» 


4.582 ±0.014 


4.27 ± 0.09 


(cgs) 


Stellar radius 


R, 


0.79 ± 0.02 


1.17 ±0.13 


Rg 


Stellar mass 


M, 


0.87 ± 0.03 


0.94 ±0.10 


M 


Semi major axis 


a 


0.0272 ± 0.0003 


0.0561 ±0.0020 


AU 


Orbital inclination 


i P 


85.35 ± 0.20 


87.4 ± 1.5 




Stellar reflex velocity 


Ki 


0.0843 ± 0.0030 


0.1101 ±0.0041 


kms s 


Orbital eccentricity 


e 


(fixed) 


(fixed) 




Planet radius 


Rp 


1.27 ± 0.03 


1.37 ±0.20 


Rj 


Planet mass 


M p 


0.46 ± 0.02 


0.89 ± 0.07 


Mj 


Planet surface gravity 


l°g£ P 


2.81 ± 0.03 


3.03 ±0.12 


(cgs) 


Planet density 


Pp 


0.22 ± 0.02 


0.34 ±0.14 


Pi 


Planetary equilibrium temperature 


T f 


1315 ± 35 


1270 ± 80 


K 


Parameter 


Symbol 


WASP-59 


WASP-60 


Units 


Transit epoch (HJD-2450000.0) 


T 


5830.95559 ± 0.00053 


5747.0295 ± 0.0022 


days 


Orbital period 


P 


7.919585 ±0.000010 


4.3050011 ±0.0000062 


days 


Transit duration 


h 


0.102 ±0.002 


0.139 ±0.005 


days 


Planet/star area ratio 


(RJRJ 2 


0.0169 ±0.0008 


0.0060 ± 0.0006 




Impact parameter 


b 


0.29 ±0.18 


0.37 ± 0.23 


R 


Scaled stellar radius 


RJa 


0.041 ± 0.003 


0.100 ±0.011 




Stellar density 


P. 


3.11 ±0.64 


0.72 ± 0.20 


Pg 


Stellar surface gravity 


logg. 


4.72 ± 0.06 


4.35 ± 0.09 


(cgs) 


Stellar radius 


R 


0.613 ± 0.044 


1.14 ± 0.13 


Rg 


Stellar mass 


M, 


0.719 ±0.035 


1.078 ±0.035 


M 


Semi major axis 


a 


0.0697 ±0.0011 


0.0531 ±0.0006 


AU 


Orbital inclination 


h 


89.27 ± 0.52 


87.9 ± 1.6 




Stellar reflex velocity 




0.1096 ±0.0044 


0.0608 ± 0.0038 


kms s -1 


Orbital eccentricity 


e 


0.100 ±0.042 


(fixed) 




Argument of periastron 


CO 


74 ± 15 






Planet radius 


Rp 


0.775 ± 0.068 


0.86 ±0.12 


Rj 


Planet mass 


M p 


0.863 ± 0.045 


0.514 ±0.034 


Mj 


Planet surface gravity 


l°g£p 


3.52 ± 0.08 


3.19 ±0.12 


(cgs) 


Planet density 


Pp 


1.8 ±0.5 


0.8 ±0.3 


Pi 


Planetary equilibrium temperature 


T f 


670 ± 35 


1320 ± 75 


K 



plotted in Fig. [7] as examples. We also estimate the stellar ages 
from the isochrones in these HR diagrams; although their uncer- 
tainties are of order several Gyr, they remain in broad agreement 
with gyrochronological ages. 

The best fits of the WASP-52, 58, 59, and 60 systems are 
plotted over the data in Figs. [1] [2] [3] and [4] The median val- 
ues and 1-cr uncertainties of the parameters obtained from the 
MCMC are reported in Table [3] Each system is detailed below. 

3.2.2. WASP-52 

There is a good agreement between the stellar mass obtained 
from the MCMC and that estimated from mass tracks in evo- 
lution diagrams for WASP-52. The red noise levels in the light 
curves are negligible; averaging the photometric data in 6- to 
10-minute bins made no appreciable difference to the error bars 
on the derived parameters. The radial velocity fit provide sys- 



temic radial velocities of -0.8953 ± 0.0075, -1.0415 + 0.0033, 
-0.9265 + 0.0062, and -0.8877 + 0.0040 km s" 1 for SOPHIE 
HE1, HE2, HR, and CORALIE, respectively, with a dispersion 
of the residuals of the order of 13 ms -1 . We fit the data using 
only parts of the radial velocity datasets, which provide similar 
results. This suggests that the systematic shift of 146 + 8 ms"' 
between the two SOPHIE HE datasets (see Sect. 12.2b has no sig- 
nificant impact on the derived parameters of the system. There 
is a hint of a drift of ~ +40 ms -1 over the 15-month time span 
of our observations. It could be due to an extra component in the 
system, but the unknown systematic shifts between datasets do 
not allow any firm conclusion to be drawn on the source of the 
drift. No significant eccentricity is detected. 

The orbital period is 1.75 day. The one-day alias of that 
signal is 2.33 days; it should be considered as the photometric 
and radial velocity observations were made with a typical 1- 
day frequency. The periodograms show signal at both frequen- 
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5000 4500 4000 6500 6000 



T eff A T eff /K 

Fig. 7. HR diagrams for WASP-52, 58, 59, and 60 (Yonsei- Yale stellar models). Each host star is denoted by a blue cross as a function 
of the stellar density and the effective temperature. Solid, black lines are the isochrones, and red, dashed lines are the evolutioanry 
tracks for different masses. WASP-52 (upper, left): 0. 1 Gyr isochrone, and then 1 Gyr to 15 Gyr in 1 Gyr increments, and mass tracks 
from 1.0 to 0.7 M in 0.1 increments from left to right. WASP-58 (upper, right): 1 to 20 Gyr in 1 Gyr increment isochrones, and 
mass tracks from 1.1 to 0.7 M in 0.1 increments from left to right. WASP-59 (lower, left): 0.1, 0.4, 0.6 Gyr isochrones, then 1 Gyr 
to 10 Gyr in 1 Gyr increments, and 0.7 and 0.6 M mass tracks from left to right. WASP-60 (lower, right): 1 Gyr to 10 Gyr in 1 Gyr 
increments, and mass tracks from 1.3 to 0.70 M Q in 0.1 increments from left to right. 



cies 1.75 and 2.33 days. However, the radial velocity data folded 
at 2.33 days clearly provide a poorer fit, with a dispersion of 
30ms _1 in the residuals, more than two times larger than the 
1.75-day solution. This confirms that 1.75 day is the correct or- 
bital period of the planet. 

WASP-52b is a hot jupiter which is less than half the mass 
of Jupiter, but with a significantly larger radius. It is thus a new 
example of an inflated hot jupiter. Together with the small radius 
of its K2V host star, it makes the planetary transits particularly 
deep at 2.7 %. The case of WASP-52b should act as a warning, 
since many transit signals of this large depth detected by photo- 
metric surveys may be flagged as due to binaries, and a priori 
removed from radial-velocity follow-up (see also Southworth 
et al. I2012I I. Hence, there may be an underlying population of 
heavily inflated planets with large transit depths that have been 
mistakenly identified as binaries and hence their true nature has 
never been revealed. 



3.2.3. WASP-58 

WASP-58 has larger uncertainties on its derived parameters as 
no extra transit light curves were obtained in complement to 
the original SuperWASP data. Assuming a main-sequence mass- 
radius prior would reduce the uncertainties, in particular on the 
stellar and planetary radii, but we did not make that assumption 
in our final results in order to remain conservative. This does not 
change the parameters significantly. 

For that target, the empirical estimation of the stellar mass 
derived from T e s and [Fe/H] (Torres et al. 120101 Enoch et 
al. 1201 II) provides M* = 1.03 + 0.05 M , whereas the isochrones 
favour lower values around M, = 0.82 ± 0.06 M . This discrep- 
ancy may be due to the low metallicity of the star, which could be 
beyond the validity of the Torres-Enoch calibration. Metal-poor 
main-sequence stars of a given mass do indeed burn hotter than 
their metal-rich counterparts. The metallicity of WASP-58 has to 
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be increased to at least -0.2 to allow a better agreement between 
the two estimations. However, when we consider other cases 
where hot jupiters orbit low-metallicity host stars (e.g. WASP- 
21, 37, or 78) we do not find similar discrepancies between the 
two estimates. Thus we remain uncertain of the true source of the 
disagreement. To reflect this, we have conservatively adopted an 
intermediate value for the stellar mass with an increased uncer- 
tainty to M, = 0.94 + 0.10M o . 

The radial velocities have a systemic value of -28.6631 ± 
0.0033kms~' and their residuals show a ll-ms~' dispersion, 
in agreement with the expected uncertainties on their measure- 
ments. No structure or drift are seen in the residuals, providing 
no evidence of any extra component in the system. 

Regardless of whether we assume or not a mass-radius re- 
lationship for the host star, we find that WASP-58b is a typical 
inflated hot jupiter, with a radius 1.4 times that of Jupiter, and a 
similar mass. Its orbital period of 5.0 days is located in the higher 
tail of the hot-jupiter period distribution. It orbits a particularly 
metal-poor star, whereas most hot jupiters preferentially orbit 
metal-rich stars. 



3.2.4. WASP-59 

The systemic radial velocity of WASP-59 is -56.7976 + 
0.0043 km s _1 with a lO-ms -1 dispersion around the circular 
Keplerian fit and no detection of any additional drift. The stel- 
lar mass estimate from the MCMC analysis agrees with that de- 
rived from evolutionary tracks in the density-temperature plane. 
As expected for late K-type stars, the isochrones are too closely- 
spaced to allow them to be distinguished, and gyrochronology 
provides the only age estimate. By allowing the eccentricity to 
vary freely, we found e = 0.10 ± 0.04 and a dispersion of the ra- 
dial velocity residuals decreased to 5 m s _I . If the marginal (2 cr 
) detection of an eccentric orbit is confirmed, then it may be in 
part due to the 7.9-day orbital period, which is relatively long 
compared to typical hot jupiters. Indeed, tidal circularization is 
less effective at larger orbital separations, and initial eccentric- 
ities in the planetary orbit may not yet have been damped (e.g. 
Hut|T981). Following Matsumura et al. (2008 ), we estimated the 
circularization time to be 5 Gyr, which is about ten times longer 
than the gyrochronological age. 

WASP-59b is one of the rare hot jupiters found by ground- 
based transit searches with orbital periods longer than seven 
days, together with WASP-8b, HAT-P-15b, and HAT-P-17b. 
Detecting such transiting objects from ground-based photome- 
try from a single site remains difficult because of the night/day 
duty cycle. The space-based missions CoRoT and Kepler pro- 
vide continuous photometry, and are therefore more suited to 
detect planets with longer orbital periods such as CoRoT-4b, 6b, 
10b, 20b, and Kepler-28c, 29b, 29c, 32c. 

The 0.7-Rj U p radius of WASP-59b places it just between the 
hot jupiters (with radii typically between 1 .0 and 1 .4 Rj up and 
periods below 10 days) and the low -radius planets recently de- 
tected mainly by Kepler (with radii typically between 0.1 and 
0.5 Rjup). Only a few objects are known in that intermediate 
radius range. One of the few others is the planet HD 149026b 
(Sato et al. 120051 1. whose density indicates a large dense core . 
Interestingly enough, the planet WASP-59b also has a particu- 
larly high density. Its host star is also one of the rare stars cooler 
than 5000 K that harbour a hot jupiter. 



3.2.5. WASP-60 

The systemic radial velocities of WASP-60 are -26.4967 + 
0.0103 km s- 1 and -26.6288 ± 0.0054km s" 1 for the two HE 
datasets. The systematic shift between them (132 + 12ms _1 ) 
agrees with that found above for WASP-52 (146 + 8 m s _1 ); this 
reinforces our conclusion that this shift has an instrumental ori- 
gin (see Sect. 12.2b . The radial velocity residuals show a 20-m s _1 
dispersion, in agreement with the expected error bars on their 
measurements. There is a hint of a drift of ~ +70 m s _I over the 
16-month time span of our observations. Again, this may sug- 
gest the presence of an additional body in the system, but the 
unknown systematic shift between the two datasets prevents any 
firm conclusion from being drawn. 

The stellar mass obtained from the MCMC analysis agrees 
with that estimated from evolutionary tracks in the density- 
temperature plane diagrams for WASP-60. The planet appears 
to have a Saturn-like radius and sub-Jupiter mass, making it an 
interestingly compact for a planet orbiting a slightly metal-poor 
star. 



3.3. Rossiter-McLaughlin anomaly analysis of WASP-52 

The radial velocities of WASP-52 measured during the 2011 
August 21 transit were fitted in order to derive the sky -projected 
angle A between the planetary orbital axis and the stellar ro- 
tation axis. We applied the analytical approach developed by 
Ohta et al. (2005) to model the form of the Rossiter-McLaughlin 
anomaly, which is described by ten parameters: the stellar limb- 
darkening linear coefficient e, the transit parameters R p /R*, a/R* 
and i p , the parameters of the circular orbit (P, To, and Ki), the 
systemic radial velocity, and finally vsin/, and A. We adopted 
e = 0.82 computed by Claret (2004) for the g' filter correspond- 
ing to the SOPHIE wavelength range. The transit and orbital pa- 
rameters were determined above from the light curves and radial 
velocity fits. Their uncertainties are negligible here for the fit 
of the Rossiter-McLaughlin anomaly. We checked that varying 
them in their 1-cr intervals has no impact on the quality of the 
fit nor on the derived A value. The main parameters that play a 
role in this fit are the systemic velocity of that particular dataset, 
A, and vsin/,. As these parameters could be correlated in the 
Rossiter-McLaughlin fit, we computed the x 2 of the fit on a 
three-dimensional grid scanning all their possible values. 

The systemic velocity was determined with a +5-ms~ 1 ac- 
curacy. It is constrained thanks to the observations secured im- 
mediately before and after the transit. The confidence interval 
contours estimated from^ 2 variations for the A and v sin z, show 
no correlation between these two parameters, as expected for 
systems with impact parameter b significantly different from 
zero. From^ 2 variations, we obtained A = 24°^ 7 and vsin/» = 
2.5 + 1 .0 km s~ l . The best fit is plotted in Fig. [8] 

The dispersion of the residuals about the fit is 20 ms -1 , 
in agreement with the expected uncertainties. The Rossiter- 
McLaughlin anomaly detection is noisy and only spans a few 
points. As a statistical test for the detection, we computed the 
X 2 over the measurements secured during the transit night for 
the fits first including the Rossiter-Mclaughlin anomaly, and then 
compared this to the fits achieved when the RM anomaly is ig- 
nored. The F-test indicates that the x 2 improvement with the 
Rossiter-McLaughlin fit has a probability > 65% to be actually 
due to the Rossiter-McLaughlin anomaly detection. We there- 
fore concluded that the spectroscopic transit is significantly de- 
tected. 
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Orbital phase 



Fig. 8. Spectroscopic observation of the transit of WASP-52b. 
Top: SOPHIE radial velocity measurements as a function of 
the orbital phase (diamonds), Keplerian fit ignoring the tran- 
sit (dotted-dashed line), and final fit including the model of the 
Rossiter-McLaughlin anomaly (dashed line). The vertical dot- 
ted lines show the times of mid-transit, first, and fourth contacts. 
Bottom: Residuals of the final fit. 



Due to the quality of the the data, we remain cautious with 
respect to the numerical result on A. The orbit of WASP-52b is 
apparently prograde and slightly misaligned. We note that the 
v sin i, value obtained from the fit agrees with the value obtained 
from spectral analysis (see Table[TJ. 

4. Conclusion 

We have presented the discovery of four new transiting exo- 
planets, WASP-52b, WASP-58b, WASP-59b, and WASP-60b, 
detected mainly from SuperWASP-North and SOPHIE observa- 
tions. As with most of the transiting planets found from ground- 
based observatories, they can be ranked in the hot-jupiter class. 
Each of them however has some peculiarities. WASP-52b and 
WASP-58b are new cases of inflated planets, increasing the pop- 
ulation of giant planets with abnormally large radii whose ori- 
gin is not well understood yet despite the numerous models pro- 
posed to explain them (see, e.g., Fortney & Nettelmann 2010). 
WASP-52 is a new case of a possibly slightly misaligned sys- 
tem; it does not fit however with the possible correlations of 
A with T eS (Winn et al. l20T0l or M p (Hebrard et al. I20TT1) . 
The hot jupiter WASP-58b orbits a particularly metal-poor star. 
WASP-59b has an orbital period longer than typical transiting 
hot jupiters, and a high density, suggesting a large core. Finally, 
WASP-60b presents shallow transits of 0.6 %, illustrating the ca- 
pabilities of ground-based surveys to detect small-size planets. 
These new planets provide new targets for follow-up studies and 
increase our statistics on well-characterized exoplanets, allow- 
ing progress to be made in the comprehension of the formation 
and evolution of planetary systems. 
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Table 4. Radial velocities (electronic table). 
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f : measurements corrected for moonlight pollution. 
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